Scaling theory of the Kondo screening cloud 

Erik S. S0rensen 
Department of Physics, Indiana University, Bloomington, IN 4^405 

Ian Affleck 

Department of Physics and Canadian Institute for Advanced Research 
University of British Columbia, Vancouver, BC, V6T IZl, Canada 

(February 1, 2008) 

Abstract 

A scaling form for the local susceptibility, derived from renormalization group 
arguments, is proposed. The scale over which the uniform part of this scal- 
ing form varies can be viewed as a definition of the Kondo "screening cloud" 
~ ^K- The proposed scaling form interpolates between Ruderman-Kittel- 
Kasuya-Yosida (RKKY) results in the high temperature limit, T ^ Tk, and 
Fermi liquid results in the low temperature, long-distance limit, T <^ Tx, 
r ^ (^K- The predicted form of the Knight shift is longer range at low tem- 
peratures where the screening cloud has formed, than at high temperatures 
where it has not. Using weak and strong coupling perturbation theory com- 
bined with large scale density matrix renormalization group (DMRG) results 
we study the validity of the finite size version of the scaling form at T = 0. 
We explicitly extract a length scale proportional to the Kondo length scale, 
^K- The numerical results are in good agreement with the proposed scaling 
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form and confirm the existence of the Kondo screening cloud. 
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I. INTRODUCTION 



The Kondo effect is probably one of the most well-studied phenomena in condensed 
matter physics. The highly successful theoretical approaches include Wilson's numerical 
renormalization group (NRG) methodS, simpler and more physical renormalization group 
(RG) approaches of Anderson et al.i and Nozieresi and exact Bethe ansatz results^. Thus it 
is perhaps surprising that a fundamental aspect of this problem remains mired in controversy. 

From the RG viewpoint, the Kondo effect is associated with very large distance scales, 

~ ae^/'''^, where a is the lattice spacing, p the electronic density of states at the Fermi 
surface and J the Kondo coupling. This scale is essentially ~ vf/Tk ~ aEp/Tx, where 
vp is the Fermi velocity, Ep the Fermi energy and Tk, the Kondo temperature, is the energy 
scale associated with the Kondo effect. Since Kondo temperatures are normally of order 
10 's of degrees, this scale is normally thousands of lattice spacings (ie. microns). A heuristic 
description of the RG results on the Kondo problem says that a cloud of electrons of this 
order of magnitude surrounds the impurity spin, forming a singlet with it. The remaining 
low-energy electronic excitations outside the screening cloud do not "feel" the impurity 
spin. Rather the screened complex acts like a potential-scatterer for these electrons, with 
a unitary limit phase shift of 7r/2 right at the Fermi energy. The largeness of this Kondo 
length scale, , in experimental systems in which the Kondo effect is apparently observed 
is rather disconcerting. Even a very dilute system with 1 part per million of impurities has a 
typical inter- impurity separation of about 100 lattice spacings, much smaller than ^k- Thus 
each impurity has many other impurities inside its screening cloud and it is surprising that 
the single-impurity Kondo effect is observed at all. Nonetheless, the impurity resistivity, 
susceptibility, etc. are observed to be linear in impurity concentration and these quantities 
seem to fit theoretical expectations. 

Comparatively little of the theoretical work on the Kondo effect has focussed on spatial 
correlations. These seem to be difficult to obtain using Wilson's method and impossible 
from the Bethe ansatz. Perturbative calculations have been performedi^ as have calcula- 
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tions using the "Nagoaka equations" a type of mean field theory!. Renormalization group 
approaches have been developed by Chen et al. § and Gani. However, Chen et al.tj only 
consider short-range correlations, with r -C and do not address the issue of the size of the 
screening cloud. We incorporate the perturbative results into our discussion of the renor- 
malization group and scaling, but come to rather different conclusions than Gani about the 
scaling variables and the size of the screening cloud. Related theoretical work has addressed 
the screening cloud in the Anderson model using various approximate methods0Hli. Some 



of the previous results have been reviewed in Ref. and \L5 . 

Experiments which have attempted to look for this large screening cloud have obtained 
mixed resultffl. In particular, the NMR experiments of Boyce and Slichtei§ were inter- 
preted to indicate the absence of this cloud. This has led to some theoretical discussion 
about the circumstances under which this cloud can be observed and even to some doubts 
about its existence. 

The purpose of this paper is to examine in more detail, the behavior of the Knight shift 
(ie. the electronic spin polarization by an applied field) in the vicinity of a magnetic impurity. 
In the next section we make a scaling hypothesis about this quantity based on standard RG 
arguments and assuming the existence of a large screening cloud. We point out the rather 
unintuitive result that the Knight shift is actually longer range at low temperatures, where 
the screening cloud has formed, than at high temperatures where it has not. We argue that 
the NMR experiments are not necessarily in contradiction with our scaling form. The basic 
problem with the experiments, according to this view, is that they only probed very short 
distances, < 3a, whereas C,k is presumably thousands of times larger than a. 

We then test our scaling hypothesis numerically. This is done using a one-dimensional 
tight-binding model. We don't expect that the reduced dimensionality is important since 
the Kondo problem is intrinsically one-dimensional anyway. Taking a spherically symmetric 
dispersion relation and a 5-function Kondo interaction, we may decompose the electronic 
degrees of freedom into spherical harmonics. Only the s-wave interacts, and this corresponds 
to a one-dimensional problem. The numerical method we use basically restricts us to T = 



and a finite length, L < 40 — 50. As will be seen, this finite length plays essentially the role of 
an inverse temperature in our scaling arguments. Both the reduced dimensionality and the 
finite length may be directly relevant to more recent experimentJii which have attempted 
to find the screening cloud using small samples with lengths of order C,k ot smaller. Recent 
theoretical work has also addressed these issue M. 

In Section |I| we present the scaling form for the local susceptibility. Renormalization 
group arguments are given bridging a high temperature RKKY form with the low temper- 



ature long distance form in a single scaling expression. Sections |T| and |^ briefiy discuss 
the form of the Hamiltonian we use in our numerical work and some details of the numerical 



method. Weak and strong coupling perturbation results are presented in section and |VII 
along with DMRG results. Finally in Section [VII1| the cross over region is studied and the 
scaling form is tested. 



II. RENORMALIZATION GROUP ARGUMENTS 

We consider the standard Kondo model: 

H = Y: eki^tAa + JSim, ■ E i^l^^A',- (2.1) 

k k,k' ^ 

Here V'ka is the annihilation operator for conduction electrons of momentum k, spin a. In 
the following we suppress spin indices which are implicitly summed over. Simp is the impurity 
spin operator of magnitude s = 1/2. The o""'s are Pauli matrices, and we set h = 1. The 
total spin operator is: 

Stot = Simp + EV^k|V'k- (2.2) 
k ^ 

(We assume, for simplicity, equal g-factors for the impurity spin and conduction electrons.) 
The Knight shift is proportional to the local susceptibility: 

x(r,T) = (l/T) < ^\t)^^{t)S^,, > . (2.3) 
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This consists of a bulk part, the usual Pauli susceptibility, p/2, where p is the density of 
states per spin, together with a local part arising from the impurity. The total change in 
the susceptibility due to the impurity, usually called the impurity susceptibility, Ximp, is: 

Ximp = / d'r[x{r) - p/2] + (1/T) < S^^^Sl, > . (2.4) 

Lowest order perturbation theory gives the RKKY result, which becomes, at rkp ^ 1: 

= f + 16r^..sth-- 2^-- (2.5) 

Vp 

Here A is the dimensionless coupling constant, 

A = pJ. (2.6) 
In the limit, r <C Vp/T, (but still r ^ l/kp) this gives the well-known RKKY expression: 

^32^'°''^^'- ^'-^^ 

A crucial feature of the Kondo problem is that, for antiferromagnetic coupling, the Kondo 
coupling increases under renormalization as the energy scale is reduced. The lowest order 
renormalization group equation: 

dA/dlnA = -A^ (2.8) 

is obtained. Here A is the momentum space cut-off, or effective band-width. This gives the 
effective coupling at momentum scale A: 

^'f'^' = l-Aln(Ao/A) - 

Here Aq ~ 1/a is the bare cut-off and A is the bare coupling (defined at that scale). The 
Kondo length scale is defined from the momentum scale at which the effective Kondo cou- 
pling constant diverges: 

a = vf/Tk ^ Ao 'eV\ (2.10) 
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A finite temperature acts as an infrared cut-off on perturbation theory so that Eq. ( p.9|) 
with A replaced by T /vp can be used to define a temperature-dependent effective couphng. 

Eq. ( |2.5| ) is only valid at high temperatures and weak Kondo coupling. As the tem- 
perature is lowered, the effective Kondo coupling increases so higher order terms become 
important. The corrections of O(A^) has been calculated. From Eq. (^.21), we see that x(r, T) 
is a sum of two terms: 

X(r) = [dr < ^t(r,0)^^(r)5Lp(r) > 

+ f dT< ^t(r)^^(r,0) / d\'^\v',T)^^{v\T) > 
Jo 2 J 2 

= Xdeir) + Xeir). (2.11) 

We have adopted the notation of Ref. ( §). Note that the sum of these correlation functions 
is independent of r since the total spin is conserved, so that the r-integral simply gives a 
factor of (3. However, the individual correlation functions depend non-trivially on r. We 
could equally well add the equal time correlation functions and multiply by P, but we choose 
the above representation because both terms have been evaluated explicitly in the literature, 
in a convenient form. The RKKY term of 0(A) in eq. ( p^.S] ), comes entirely from Xde- In the 
asymptotic region, rkp ^ 1, Ep/T ^ 1, with r <^ Vp/T, including the correction of O(A^): 

X{r) -P- = + X'MkFT) + constant]}. (2.12) 

The logarithmic term comes entirely from Xde- It was first calculated in Ref. ( ^ and ( 
Xe contributes only to the constant in Eq. ( p.l2|) . [See Eq. (B2) of Ref. ( H).] 



Note, from Eq. (|2.9| ), that the quantity in brackets in Eq. ( |2.12| ) may be written AefT(r) + 
AefT('")^ + constant, to O(A^). This expression exhibits an infrared divergence at large r. That 
is, for sufficiently large r, r > C,k, the O(A^) term exceeds the 0(A) term. Note however, 
that this correction term is at least finite as T — 0. ie. it is Acfr(r) that appears, not 
XesiT). Thus, at least to O(A^), a finite r is acting like a cut-off on the infrared divergences 
of perturbation theory. It is an important question whether or not this persists to higher 
orders in perturbation theory, ie., is perturbation theory valid for r <^ C,k even for T <^ T^, 



with the actual expansion parameter being AcfT(r)? Based on an examination of higher order 
terms GanI has argued this not to be the case. He claims that higher order terms diverge 
as T ^ for non-zero r and that it is therefore necessary to have T ^ Tk for perturbation 
theory to be valid. [This point will be examined in detail in Ref. ( ^0|).] In this case, it is 
probably more useful to rewrite Eq. ( ^.121 ) in terms of Aefr(T). To O(A^): 



'2 = S^^^^^^^^) + XeffimnirT/vp) + constant]}, (2.13) 
using Ep/kp ~ vp- 

At very low temperatures and large distances, T -C Tk, r ^ ^k, we expect x(^) to 
be determined by the zero-energy fixed point. Within the local Fermi liquid theoryi of 
this fixed point we can then estimate The zero-energy fixed point corresponds to 

a screened impurity which just acts as a potential scatterer for the low energy electronic 
degrees of freedom, with a phase shift of 7r/2 at the Fermi energy. The local susceptibility of 
a potential scatterer follows directly from the formula for Friedel oscillations in the electron 
density, n(r), with an s-wave scatterer and a 7r/2 phase shift. For kpr ^ 1: 

n(r) = no - cos[2/cFr + 7r/2]. (2.14) 

27r^r'^ 

Noting that a magnetic field, H, simply shifts the chemical potential by ±gfiBH/2 for spin 
down or spin up electrons, we obtain: 

Note that x(r, T) is longer-range at low T after the screening cloud has formed, x 1/r^, 
than at higher T before it has formed, x 1/t^, (Eq. ( |2.5| )). An analogous result occurs in 
spin chain systems^. 

Corrections to Eq. ( p.l5| ) can be derived by doing perturbation theory in the leading 
irrelevant operator. Part of the leading correction can be obtained by considering a field- 
dependence of the phase shift:! 

5" = 7i/2 + ahc/TK, (2.16) 



where c is a dimensionless constant of 0(1). Generalizing the Friedel oscillation formula of 
Eq. ( p.l4| ), the local density of spin-cx electrons becomes: 



n.(r) = Y - ^ cos[2kUh)r + 6^{h)]. (2.17) 

Upon differentiating with respect to h to obtain the local susceptibility, we now obtain an 
additional term: 

Sx= , cos(2A:^r). (2.18) 

Note that this term drops off more rapidly with r than the term in Eq. (|2.15|) obtained from 
differentiating kp{h) and is smaller for r ^ C,k- A very similar Fermi liquid calculation of 
xi'f'), in the Anderson model, was performed in Ref. |13[ However, this calculation effectively 
ignored the field dependence of kp, and hence only obtained the subdominant term of Eq. 
( p.l8| ), not the leading term of Eq. ( p.l5| ). Explicitly, in Eq. (2.10) of Ref. |13|, the free 



electron Green's function, Fr{iuji) must be evaluated in a finite magnetic field. Taking this 
into account, we obtain our expression, Eq. ( |2.15| ). 

According to Fermi liquid theory, the impurity susceptibility, Ximp, defined in Eq. ( p.4| ) 
is 0{1/Tk). This appears to arise from a short-range part of %('") which does not oscillate 
at wave-vector 2kF- However, since the impurity has been "integrated out" to obtain the 
Fermi liquid theory, it is difficult to ascertain how much of Ximp comes from x(r) and how 
much comes from the impurity self-correlation function, in Eq. ( p.4|) . It is also difficult to 
tell whether the contribution, if any, from x(r) has a range of 0{^k) or only of O^l/kp), 
because the cut-off has been reduced to 0(1/^a') to obtain the Fermi liquid theory. 

We now wish to formulate a scaling hypothesis for x(r, T) which we expect to be valid 
at arbitrary r and T in the scaling region, r ^ a, T ^ Ep- For this purpose it is very 
convenient to use the relativistic one-dimensional formulation of the Kondo problem. [Ref. 



22| The mapping to one-dimension is exact for pure s-wave scattering. The use of a reduced 
bandwidth and linear dispersion relation which leads to the relativistic model is expected to 
be valid in the scaling region. The three dimensional electron field is expanded in spherical 
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harmonics and then the s-wave part is written in terms of left and right moving components 
(ie. incoming and outgoing): 

V'(x) = — ^ — fe"*''^''V'L(r) - e^''^''V'ii(r)l + higher harmonics. (2.19) 
2V27rr '- 

The left and right-moving fields, defined on r > obey the boundary condition: 

MQ)-MQ)- (2-20) 

We may flip the right-moving fleld to the negative axis, so that we work with left-movers 
only deflned on the entire real axis: 

^jJJ^{-x)=^|Jn{x). (2.21) 

The one-dimensional Hamiltonian can be written: 

dri;l{r){id/dr)Mr) + VFXtl;i{0)-MQ) ' Simp- (2.22) 

-oo ^ 

X — p/2 can be expanded in spherical harmonics; only the s-wave harmonic is non-zero. This 
can be written in terms of one-dimensional uniform and 2kF susceptibilities: 

^ - = ^I^- + (e'''"''X2;t, + cc.)], (2.23) 
where cc. denotes complex conjugate and: 

Xunir.T) = (1/T) < [^l:{{r)^^|;^{r)-r^Pi{-r)^M-rWT > 
X2kAr,T) = (l/T) < V^t (r)yV'4-r)5| > . (2.24) 

Here St is the total spin in the one-dimensional theory: 

St = Si„,p + ^ r dT^{{T)^Mr)- (2.25) 

X2kF can be shown to be real using particle-hole symmetry. This follows since under 
particle-hole symmetry: 

Mr) ^ <T^V'l(r), (2.26) 
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and hence: 

— > Sj' 

iljl{r)cTiljL{-r) ^pLi-r)f^i^L{r) = tpi{r)(TiljL{-r) \ (2.27) 

When particle-hole symmetry is broken, as it is for a realistic Hamiltonian, we expect X2kF 
to have a phase, 6, which is non-zero but constant in the scaling region. This can be seen 
from spin-charge separation in the one-dimensional formulation of the Kondo problem. The 
Kondo interaction, which produces the non-trivial scaling behavior, occurs entirely in the 
spin sector. In the absence of particle-hole symmetry there is a marginal potential scattering 
term, —9ip'l{0)ipL{0), which is a pure charge operator. Upon bosonizing, this is linear in the 
charge boson and hence doesn't renormalize (is exactly marginal). We can write X2kp as 
a product of spin and charge correlation functions. The charge correlation function just 
contributes a constant factor e*^ to Xkp- 

In the lowest two orders of perturbation theory, discussed above, the function Xun vanishes 
at r ^ l/kp. In fact, it is possible to prove that this happens to all orders in perturbation 
theory.0 This is also consistent with Fermi liquid theory, given the uncertainties in that 
theory, discussed above, about the origin of Ximp- 

We expect the one-dimensional local susceptibility to obey scaling in the following sense. 

After extracting a factor of l/vp, X2kp could, in principle depend on three dimensionless 

variables, which can we taken to be, rT/vp, A and D/T. Here D is the effective bandwidth 

in the one-dimensional theory, a quantity of 0{Ef). The scaling hypothesis asserts that the 

bare coupling constant, A and D/T do not appear independently but only in the combination 

making up the renormalized coupling constant, XesiT). The dependence on XesiT) may be 

exchanged for a dependence on T/Tk- To see that these two quantities are related, note, 

from Eqs. (|2.9| ), and (|2.10| ), valid at weak X^s, 

1 

.Aeff(T)_ 

In the intermediate to strong coupling region, the value of Aefr is non-universal, ie. am- 
biguous. In this regime it is better to use T/T^ as a measure of the dimensionless effective 
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exp 



T/Tk. (2.28) 



coupling. This is in accord with the idea that I/Tk is the couphng constant for the lead- 
ing irrelevant operator at the low-temperature fixed point. As usual, we multiply by the 
effective cut-off, T, to form the dimensionless coupling constant. Thus we write the scaling 
hypothesis as: 

X2k, = —f{rT/vF,T/TK), (2.29) 

Vp 



where / is a real universal scaling function. Using Eq. ( p.23|) the equivalent statement for 
the three-dimensional susceptibihty is: 

X - p/2 = I ^ cos{2kFr)f{rT/vF, T/Tk). (r > l/kp] T, Tk < Ep) (2.30) 

oTT VpT 

As we will see below, this is consistent with what is known about the behavior at T ^ Tk 
from perturbation theory and T <^ and r ^ S,k from the local Fermi liquid description 
of the critical point. Note that this scaling hypothesis does imply the existence of the large 
screening cloud, since if T < T^, the length scale over which x(r) varies is at least 
(apart from the 2kp oscillations and the short-range part). In fact, this scaling hypothesis 
is perhaps the best definition of what it means to have a screening cloud. Note that this 
scaling form does not include any anomalous dimension. We expect this to be absent since 
Stot is conserved. The more general case, with unequal gyromagnetic ratios for electrons 
and impurity, involves a non-conserved operator. This will be discussed in Ref. ( |20|) . 

The known results from perturbation theory and Fermi liquid theory, Eq. (|2.5| ), ( p.l3D , 



( P.15| ) and ( p.l8| ), are all consistent with this scaling hypothesis and imply certain limiting 



forms for the scaling function. Eq. ( |2.5|) implies that 



/(a:,y)- 17-777-^^' (2.31) 
2 smh(27rx) my 



for y ^ 1. Eq. ( p.l3| ) gives a higher order correction in l/lny to f, when x is also small. 



Eq. (p.l5|) and ( p.lSj) imply that, for y <^ 1 and x/y > 1, 



/(x, y)^2 + constant ■ -. (2.32) 
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The function / in the regime y < 1, x/y = r /^k < 1 is of special interest. It describes 
the interior of the screening cloud at low T. One might naively suppose that a small r <^ 
would also cut off the renormalization of the effective coupling so that deep inside the 
screening cloud we recover weak-coupling behavior (for weak bare coupling) even at low T. 
As mentioned above, Gani has argued, based on higher order perturbative calculations, that 
this is not the case. It follows that the scaling function would be non-trivial in this region. 

Now let us consider the experiments of Boyce and Slichter on Fe doped Cu. They 
measured what they interpreted as the Knight shift from 5 different shells of Cu atoms at 
distances up to 5th nearest neighbor. We note that for the Cu fee lattice, assuming the Fe 
impurities occupy Cu lattice sites, the fifth nearest neighbor is at a distance of \/Qa ^ 2.4a 
where a is the nearest neighbor separation. The measurements were taken from T=300K 
down to well below what is believed to be the Kondo temperature of 29K. They found the 
factorized form: 

x(r,T) = -^, (2.33) 
1 + Ik 

for some rapidly varying function /(r) (which, in fact, changes sign over the small range of r 
considered). Note that all measurements are taken in the regime rTjvp ^ 1, rT^jvp <^ 1. 
In fact the values of r are so small that it is unclear whether the scaling form of Eqs. (p.30|) 



holds at all. In particular, the short-range part of x may be contributing. If we assume r 
is large enough that this can be ignored, and the scaling form holds, then we may consider 
the short distance limit r ^ vf/T,vf/Tk of the scaling function. According to Gani the 
behavior of f{x,y) is non-trivial in the lower temperature range of the experiment, ?/ <^ 1, 
x/y <^ 1. ie., we do not know the behavior of the scaling function at low T deep inside 
the screening cloud. It is possible that f{x,y) exhibits an approximately factorized form for 
x/y <ti 1 and all y: 

f{x,y)^^^V (2.34) 
y + 1 

This behavior would explain the experimental results. Note that such factorization could not 
also occur at large r, r ^ C,k if our assumed scaling and asymptotic behaviors are correct. In 
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this region, the behavior at T ^ Tk crosses over to at T ^ Tk- More experiments 
at larger r could clarify the situation. Experiments in the region 1/kp <^ r <^ C,k niay be 
feasible . These would probe the short-distance part of the scaling function. A full study of 
the scaling function would require going out to values of r > One point to bear in mind 
is that, assuming the existence of a large screening cloud, the average impurity separation 
is much less than C,k so inter-impurity interactions may be playing a large role. This may 
make experimental observation of the screening cloud very difficult, at least until a better 
understanding of the effect of inter-impurity interactions is obtained. 

III. TIGHT BINDING MODEL 

In order to apply the density matrix renormalization group (DMRG) method we rewrite 
Eq. in real space. This is the standard s — d Kondo model. The model is described by 
a tight binding Hamiltonian coupled to a s = 1/2 impurity spin. Simp. 

L-l 

H=-tY. (V-J^.+La + 4lli^^,c) + Hk. (3.1) 
1=1 

Here describes the coupling to the impurity spin. For the bulk of our results we consider 
a single s = 1/2 impurity spin at the left end of an otherwise open chain. In this case the 
coupling to the impurity spin described by Hk takes the form 

= \ j ((a\V^it - ^lii'n) Simp + ^It^U'^imp + V^u^iT-^ilp) • (3-2) 

We also briefly consider the case of two s = 1/2 impurities located at either end of the chain, 
correspondingly Hk becomes: 

Hk = JSi • ^< + JSl ■ A^^i^i (3.3) 
In all our results below we have used t = 1. 
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IV. NUMERICAL METHOD 



We use the density matrix renormalization group (DMRG) method as developed by 
S. R. White and R. M. NoackH. For a detailed explanation of the method we refer the 



reader to Ref. 23. The method is centered around calculating the density matrix and its 
corresponding eigenvalues and eigenvectors for two appropriately defined parts of the total 
system. Usually this is done by simply considering the system as having a left and a right 
part. The eigenvalues of the density matrices can then be interpreted as the probability for 
the subsystem to be in an eigenstate of the density matrix given the constraint that the total 
system is in a fixed "pure" eigenstate (usually the ground-state). For stable fixed points 
of the DMRG method it can be shownS that states in the thermodynamic limit are well 
represented by "matrix product ground states"iHH^. 

In an iteration the density matrices for each half of the system contain 4 x m states 
(including the S = 1/2 impurity spin). Of these states m are kept to start the next iteration. 
We have used m in the range 128-200. It is also extremely useful to use all the symmetries 
of the original Hamiltonian and constrain the whole calculation to a subspace defined by 
suitable quantum numbers. This increases the precision of the method dramatically. For 
the Kondo Hamiltonian that we consider here we have in addition to the total z-component 
of the spin, (including both the impurity and electron spin), and the number of electrons 
also parity, P, for reflection around the midpoint of the chain and, RPH a particle hole 
symmetry combined with a rotation that changes the sign of the z-component of the spin. 
The first two are diagonal in the usual product basis and are thus trivial to implement, the 
P and RPH symmetries are non diagonal and a considerable effort has to be expended to 
implement these symmetries. The parity P, is standard and takes the two values 1,-1. 
This is only applicable when we consider two impurities since the one impurity model is not 
symmetric with respect to a reflection around the middle of the chain. The RPH symmetry 
is an on-site symmetry and can be used both for the one and two impurity model that we 
consider. It is defined by 
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RPH: i,,^-. (_iy^t/3(a-)^^ 
RPH: ^f-. {-ly^a^y^^^p 

RPH : |0 > ^ |F > . (4.1) 

Here |0 > is the empty state and \F > the completely filled state. This is an exact symmetry 
of the Hamiltonian and it commutes with P. We can then specify a state by the four quantum 
numbers, the filling factor, S^, P, and RPH. In the following we shall always work at half- 
filling, but we chose the remaining quantum numbers so as to select appropriate states. 
For calculations performed on a system with an impurity at both ends we can use all four 
quantum numbers, for the case with only one impurity the chain has no longer refiection 
symmetry around the midpoint of the chain and P is no longer a good quantum number. 

The bulk of our results are obtained for the case of only one impurity. In order to 
obtain optimal precision it is necessary to use a combination of the so called "infinite" and 
"finite" length DMRG methods!!. At each step of the infinite chain method a complete 
run of the finite length method is performed. This is done in the following way. We start 
by considering a 2 site system in addition to the impurity spin. Using the infinite lattice 
method we generate a 4 site system with matrices representing the impurity spin plus the 
2 left sites of the chain and another matrix representing the two right sites. We denote this 
by 2l + 2{i, where the impurity spin is included in the matrix 2^. The finite system method 
is then used to arrive at a system consisting of 1^ + 3^? where the exact 8x8 matrix is used 
for the impurity spin and first site for the matrix 1^. At this step the expectation values 
< Sj^ip > and < Si > are calculated since these matrices are known exactly. The finite 
lattice method is then used to generate 2l + 2^ and 3l + 1r, at each step calculating the 
expectation value of the electron spins at the sites where the matrices are known exactly, 
i.e. at sites 2,3 and 3,4 respectively. Then we use the finite system method to generate 
3l + 3/j and the whole procedure is repeated. For large chain lengths this is exceedingly 
slow, but we have been able to treat chain lengths of up to 50 sites. This method is exact 
out to the point where the matrices A''^ and A^r have to be truncated and it has the great 
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merit of yielding the same precision for the expectation value of the spin components in the 
middle and at the end of the chain. This latter point is essential to our analysis of how the 
impurity spin is screened at each chain length. 

In all cases below we shall take t — 1 and we always take L even. Furthermore, we shall 
usually work in the ground-state subspace which for 1 impurity is given by the quantum 
numbers, = 1/2, RPH = 1. For 2 impurities the ground-state subspace is defined by, 

= 0, RPH = 1, and P = -1. 



V. FREE CHAIN 

Let us first consider the tight binding model in the absence of any impurity. 

L-l 



i=l 

Here the subscript "free" denotes the free chain. Since we shall be concerned mainly with 
free boundary conditions for the chain we introduce two "phantom sites", and L + 1, where 
we require that ^jJ{0) = 'ijj{L -|- 1) = 0. This model can be solved by transforming to Fourier 
space and we obtain, in units of the lattice spacing, a, 

Hfree = -2t ^ cos(aA;)^l"V'fe,„ (5.2) 

k 

where the allowed values for the wave vector, k, are 

7VTI 

k = j^, n = l...L, (5.3) 

and thus < A; < tt. We have now essentially two decoupled Fermi seas for up and down 
spin electrons respectively. The energy is then given by the expression 



-E'free — — 2t 



cos(aA;^) -|- ^ cos{ak^) 



(5.4) 



Here kp,kp denotes the Fermi wave vector for up and down spin electrons. We can now 
choose a simple representation for the field operators. 
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V'i = y^^IIsin(A;j>fc, V"] = y^^I] sin(A;j)4. (5.5) 
We note that with this definition the field operators obey the commutation relations 

This solution leads to the interesting fact that for a chain with an odd number of sites 
at exactly half filling, i.e. one electron per site, the magnetization per site (the expectation 
value of the z-component of the electron spin, < Sj >, is given by the following expression 

<s;> = l< v^jV] - V^j Vj > 

= sm — . (5.DJ 

L+l 2 ^ ^ 

Thus the on-site magnetization, < Sj >, is non-zero only on odd sites, an artifact due to 
the open boundary conditions. This result can be reproduced by the DMRG method. 

For an even length chain, at half-filling, the on-site magnetization is always zero in the 
ground-state where St = 0, since the expectation value of the electron spin is zero at every 
site due to rotation symmetry. However an excited state with 5*^ = 1 will have unpaired 
spin up electrons in the states with k = L7t/{2L + 2), (L + 2)7r/(2L + 2). Redoing the above 
calculation above for an even length chain in the St = 1 state we obtain: 



^ L + l 



sm — + sm 



2(L+1) 2(L+1) 
Equivalent results can be obtained for higher excited states. 



(5.7) 



VI. WEAK COUPLING PERTURBATION 

A. < 5| > 

By considering the term in Eq. ( p.l[ ) as a perturbation we can do first order pertur- 
bation theory in (J/t). We start with a system at half filling with an even number of sites 
L and we shall take the total z-component of the spin to be 1/2, S^ = 1/2. The impurity 
spin is at the far end of the chain at site 1. Here we consider only one impurity described 
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by Hk in Eq. ( p.2|) . Since we are only considering first order perturbation theory only the 
term involving in Hk will contribute. Thus, we have for the unperturbed ground-state 

00 = 10 >free X | t> • (6.1) 

Here |0 >free is the half-filled sea described in the preceding section and | |> denotes the 
spin up state of the impurity spin. The first order perturbation to the wave function is 

1 -P J 



>-- 



- 4^^i) 10 >free X | T> ■ (6.2) 



Eq — i/frce 4 

Here P is the projection operator onto the ground-state Eq. (|6.1| ). When L is even all single 
particle states with energy Ek = — 2t cosfc below the Fermi energy at ep = 0,kp = 7r/2 are 
filled for the unperturbed chain. With this notation we can rewrite Eq. ( |6.2|) as 



I , J v-^ sin/csinfc' / I i \ 

101 >= or I O ( «fc'T«fcT - ) 10 >free X | t> • (6.3) 

+ ^ k,k',e,<0,e,,>0 ^k-Sk' 

Thus, the on-site magnetization becomes 

2J sin fc? sin fc'? sin sin fc' , 

< S- >= 77 y . (6.4) 

If we now set k' = tc — k" we see that Sk' = —ek",sink'j = — (— !)■' sin A;"j and we can 
therefore rewrite the above equation as 

< gz J (-1)^ sin kn sin km sin kj sin k^j . . 

^ ~t{L + l)\^^^ cosfc„ + cos/c™ ' 

where as above we have kn = nn / {L + 1) . For j ^ 1 Eq. ( |6.5| ) can be analytically evaluated 
as: 

This formula is basically the ID version of the RKKY formula, with replaced by 1/r 
for trivial dimensional reasons. Note that this expression only has a staggered part, not a 
uniform part. The uniform part is O(J^). It is quite easy to see that the sum over all j of 
Eq. ( |6.5D vanishes exactly. (This is simply a consequence of the fact that the total electron 
spin of the unperturbed ground-state, with L even, is zero.) 
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The above result, Eq. ( |6.5| ), can be compared to DMRG results obtained for weak cou- 
plings. In Fig. |1] we show results for a 30 site system with one impurity at the left end. The 
circles denote < 5*1 > calculated with the DMRG method, while the crosses are Eq. (|6.5D. 
The dotted line indicates the asymptotic form Eq. ( |6.6| ). The calculation was performed 
keeping m = 128 states in the ground-state subspace defined by RPH = 1, = 1/2. 
Clearly there is a very good agreement between the perturbative results, Eq. ( |6^ , and the 
DMRG results. 

Equation ( |6.5|) can also be compared to results obtained using the DMRG method with 
two impurities, one at either end of the chain. In Fig. ^ we show the expectation value of 
the z-component of the electron spin as a function of chain index, j, for a Kondo coupling, 
J, of 0.05. In this calculation we have a S = 1/2 impurity spin at each end of the chain. 
The total z-component of the electron and impurity spin was chosen to he = 1 so as to 
polarize the two impurity spins as much as possible, m = 150 states were kept and we fixed 
the two remaining quantum numbers RPH = 1, P = 1. As usual we work at half filling. 
The DMRG results are shown as the circles in Fig. |^. In order to compare these results to 
Eq. ( |6.5D we must sum the contribution from the impurity at both ends of the chain. If we 
denote by S^^^^ the result of Eq. ( |6.5| ) for the contribution for one impurity spin we obtain 

Sf^ = Sf^ + Sf%. (6.7) 

This expression is shown as the x in Fig. 0. As clearly seen in Fig. ^ there is an excellent 
agreement between the perturbation results and the DMRG results for the chain length, 
L = 60, considered. 



The expectation value of S-^^^ can also be evaluated in first order perturbation theory. 
Using Eq. ( |6.2[ ) and remembering to include the contributions from the x and y parts of 
Kondo interaction, we obtain: 
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< >^ 1/2 



J 



2 



E 

efc<0,e^/>0 



sin k sin k' 



2 



(6.8) 



Where we must be careful to include the correction to the ground-state wave-function nor- 
malization. In the large L limit, replacing the sum by an integral, the integral has a log 
divergence. Thus the last term goes like InL. 

In Fig. ^ we show DMRG results for the z-component of the impurity spin as a function 
of chain length, L. The results shown are for a very weak Kondo coupling, J = 0.05. One 
impurity is present at the left end of the chain. The crosses indicate the first order pertur- 
bation result, Eq. (|6.8| ), the circles denote DMRG results for the state = 1/2, RPH = 1. 
A good agreement between the perturbative and numerical results is evident. 

C. Correlation Function 

It is also straightforward to obtain the correlation function < S^-^^Sj > to first order in 
J/t. It is simply 1/2 times Eq. In Fig. | we show DMRG results for < S-^^^S^ > 

(circles) for a 50 site chain with impurities at both ends. The Kondo coupling is very weak 
J = 0.05. Note that only half the chain is shown. The calculation has been performed with 
m = 150 states in the subspace RPH = 1, P = 1, and = 1. The crosses denote the 
results from Eq. ( |6.5|) multiplied by 1/2. We see an excellent agreement between the DMRG 
and the perturbative results. 

VII. STRONG COUPLING 
A. Energy 

In the J = oo limit the impurity traps an electron and forms a tightly bound singlet 
leaving a free chain with L — 1 sites. For finite J the singlet can be polarized and we can do 
perturbation theory in the hopping term between the first and second site. We thus take 
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H = JSi^p ■ ^l"^V^f + //frec(^ - 1) + ^, (7.1) 
V=-t (^I>2,a + V'^V'l.a) • (7.2) 

Here -fffree is the free chain Hamiltonian Eq. ( |5.1| ) for L — 1 sites. In the following we shall 
regard the hopping term, V", between site 1 and 2 as the perturbation. 

The strongly bound singlet on the first site of the chain can be excited into any of the 
seven excited states shown in Fig. ^. The excited states form a quadruplet with energy 
3 J/4 and a triplet with energy J relative to the ground-state. The perturbation, V, has 
only non-zero matrix elements between the ground-state and the quadruplet. Calculating 
the partial matrix elements in the impurity part of the Hilbert space we find 



^10 


=< 1\V\0 >= 


^4' 






=< 2\V\0 >= 








=< 3\V\0 >= 








= < 4:\V\0 > = 




(7.3) 



Following standard second order perturbation theory we can now calculate the energy shift 
due to the perturbation V. If we denote by |0 > \F > the unperturbed ground-state 
composed of the singlet and the free chain we find 

AE = y < F| < 0\V\a >< alVlO > \F > 

o J ^ I II II I 

■^'^ a=l 



-— <F|{^/.+",^2j|i^> 
bJ 

3J' 



(7.4) 



We thus find for the total energy for 1 and 2 impurities, respectively, 

-^1 impurity 4 -'-'freely -'^ ^) 

6 

E2 impurity = + Eiree{L — 2) — . (7.5) 

These perturbation results compare favorably with the DMRG results at sufficiently strong 
coupling. 
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B. wave- function 



We can now calculate the wave-function to second order in perturbation theory. We find 

\^ >= |o > |F > -— ^ Ko|« > > +77(7)' E ^5aV;o|5 > |F > + ■ ■ ■ . (7.6) 

For the calculations we shall consider here these are the only terms that will contribute. We 
need to calculate V^s. We find for the partial matrix elements in the impurity part of the 
Hilbert space, 



Vis 


=< i\V\h >= 








=< 2\V\5 >= 






V^35 


=< 3|V^|5 >= 


-^^'^ 




V^45 


=< 4|V|5 >= 




(7.7) 



We can now proceed to evaluate the expectation value of S^, 5*1, S^^^ in perturbation theory. 

f-~l QZ CZ 

'-'j ' '-'imp 

We begin by considering the z-component of the electron spin on the first site, 5*^^. In 
the unperturbed system must be zero since it is locked in a singlet with the impurity 
spin. The first non-zero contribution to Sf is second order in t/J. Since 

S',\0>=-^\5>, (7.8) 

we find that only the second term in Eq. ( |7.6| ) will contribute and we get 

< ) >= 2 < F| < 01^,^(4)^ E > \F >= -U-y < F\4'4 - 4'4\F >(^) . 

(7.9) 

From this we can derive two results; if the state \F > describes an unperturbed chain with 
L — 1 sites we find from Eq. ( ^.6|) 
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< Sf' >= -lijf^, (7.10) 

since in that case the matrix element in Eq. (|7.9|) is simply 2/L. We can also consider the 
case where \F > describes a free chain with L — 2 sites but now in a state with = 1. This 
is convenient for comparing with results with 2 impurities where 2 sites will be quenched 
out. It is in that case convenient to work in the state with = 1. We then find 

< Sf^ >= -^l)^^[sin2 4^ + sin^ V^]- (7-11) 
^ 3 L-V 2(L- 1) 2(L- 1)^ ^ ^ 



In a similar fashion can be calculated to second order int/J. We first consider the 



.2(2) 

imp 

case with one impurity. Again S'f is zero in the unperturbed state. For we find 



imp ^-^'^ ^ .^^^ ^ "imp 

5imp|0>=^|5>. (7.12) 

Thus we get the same term as before but with a different sign. However, since S'j^p is non- 
zero in the quadruplet we get an additional term from the second term in Eq. ( [7.6|) . We 
then have 



a=l 



= ^-^{jY<F\4'4-4'4\F>^'^ . (7.13) 

As before we can now obtain the results for the case where \F > describes a chain with L — 1 
sites, ie when only one impurity is present. 



And equivalently for the state with L — 2 sites and = 1 where we have the 2 impurity 
case in mind. 



<^S>-^5(« ,,15, 



To this order in perturbation theory we have in addition the following equality for the second 
order contribution to < 5*1 >: 

24 



< sf^ >=-< Sf^ >-< > . (7.16) 

Note that for < > we still have the term 1/L from Eq. (|5.6| ) which we can apply here at 
strong coupling. These equations are obeyed to a high accuracy at sufficiently high coupling 
J. To illustrate this we show in Fig. ^ results for one impurity with a Kondo coupling of 
J = 10 for L = 24. The calculation is done keeping m = 128 states with RPH = 1 and 
= 1/2. The circles denote the DMRG results for < > and the square denote < S^^^ >. 
The crosses denote the results from Eq. ( p.6| ) for a 23 site chain with the corrections from 
Eqs. ( [7.11| ) and ( |7.16| ). The plus is the result for < S'j^p > from Eq. ( [7.14] ). From the results 



in Fig. P we find excellent agreement between the above perturbation results and our DMRG 
results at strong coupling. 

We note that the fact that < Sj constant/L for j ^ ^k, and thus in a sense is 

longer range than RKKY (1/r), is analogous to the prediction that xi'^^j T) should be longer 
range when the screening cloud has formed (see discussion below Eq. (|2.15|) ), than when it 
has not. 

D. correlation function 

It is also rather straight forward to calculate the correlation function < S-^^^^Sj > to 
second order in t/J. Again the unperturbed result is just zero for j > 1. Let us first 
consider the case where j > 2, in which case we get from Eq. ( [7.13| ) 



on + 

= y(jmL-2,j-l) j>2. (7.17) 

Here we consider the case where \F > describes a free chain with L — 2 sites corresponding 
to two impurities, and G is given by 

G(L, j) =< FK^IVI - ^lV{)(^jVj - i^M)\F > . (7.18) 

For the states describing the free chain G can easily be calculated. For the state that we 
consider here with L — 2 sites and S"!^ = 0, we find 
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•-^(^'•^^ " rnru2^ sinA;sinfcj)( ^ sin A; sin /cj) 

v-^ + ^) k<kp k>kp 



( ^ sinA;sinA;j)l (7.19) 



As before we have /c = 7™/(L + 1), ri = 1 ■ ■ ■ L. When j = 2 we have to be somewhat more 
carefuL Reanalyzing the two terms in Eq. ( |7.13| ) we find that due to the non-commutativity 
of S"! and V we effectively get a sign change on the first term and thus 

< SLp^l >(^^= l{^fG{L - 2, 1) = \{^f. (7.20) 

Furthermore we must have that 

< S?Sl >(2)= _i _ (7.21) 

imp 1 4 3(/ 



for the case where we have only one impurity where AE is given by Eq. (|7.4| ). Except for 
the first argument to G we note that the above results for the correlation function do not 
depend in any essential way on whether we consider 1 or 2 impurities. 

We illustrate this with a calculation of < S^^^S^j > for a 30 site chain with an impurity 
at both ends. We keep m = 150 states and work in the ground-state subspace with RHP = 
1,P = -1, and = 0. We first consider < S^^pSf >. From the DMRG we find < 
'5'imp'S'i >= —0.24636. This can be compared with the result from Eq. ( |7.21| ) which gives 
< 5f^p5i" >(2)= -.24555 in good agreement. The results for j > 2 are shown in Fig. [^. The 
circles denote the DMRG results and the crosses the results from Eqs. (|7.19|) and ( [7.20 ). 



Clearly the discrepancy is largest for j = 2, 3. As a function of chain length, L, we have 
observed that this discrepancy decreases. Over all the agreement with the perturbative 
results is very good. 

VIII. CROSSOVER REGIME 

Having checked that the DMRG method yields the expected perturbation results in the 
strong and weak coupling limit we now proceed to study the crossover behavior and the 
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scaling predicted by Eqs. ( |2.29|) . Since the numerical results are all obtained at T = for 
finite systems, where S^^^ is a constant, the role played by Ax = x(r, T, J) — p/2 is taken 
by the the on-site magnetization, AS =< Sj > (expectation value of the z-component of 
the electron spin at site j). Hence, a generalized finite size form of Eqs. ( |2.29| ) should apply 
to < 5*1 > at T = 0. Such a finite-size form, applicable to the numerical results for finite 
L, is easily obtained; we simply substitute L for the thermal length, vp/T. In this way we 
get rT/vp ^ r/L and T/Tk — > ^k/L. We note that in order to apply Eq. (|2.29|) which was 
derived for a susceptibility (Eq. ( p.3|) ), to the on-site magnetization we need to multiply by 
T. We then obtain: 



where we have kept the finite-size equivalent of Xun as well as X2kp- We note that this form 
assures that < Sj >~ {C,k/L) J dxg{x, L/^k) is a function of C,k/L only. 

A. Scaling of L < 5| > 

We now proceed to test the above scaling form for the finite systems we have been able 
to study numerically. We start by considering how the weak and strong coupling results 
of Sections and [VI 1] can be cast into a from consistent with Eq. ( ^.1|) . First we note 
that the strong coupling expression Eq. (|5.6|) obviously obeys the scaling form. However, 
the asymptotic weak coupling result, Eq. ( |6.6|) doesn't seem to obey the scaling form due 
to the explicit dependence on J/t. Fortunately it is possible to remedy this by noting that 
J/t ^ l/ln(^^/L). This we see in the following way: We can write our 2 scaling variables 
as j/L and ^k/L. Alternatively, and perhaps better, we can replace ^k/L by the effective 
renormalized dimensionless Kondo coupling at scale L. Here we define a dimensionless 
Kondo coupling as A = J/t. Combining Eqs. ( |2.9| ) and ( |2.10| ) we get using A^^ = L: 
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As L becomes small (ie. approaches 1) \eff{L) approaches the bare coupling constant, A. 
Thus, for small L, and weak coupling, we have, by substituting l/ln(,^i^/L) for {J/t): 

perfectly consistent with scaling. 

First we consider L < S'£/2 >• this case the scaling function should be 

g{l/2,L/^K{J)) + {—^)'"^'^f{^/'^iL/^K{J))-, and we can suppress the first argument. In 
addition since we only consider the case where L/2 is even, the scaling relation takes the 
simpler form L < 5*1/2 >= ^{L/C^k)- In Fig- H we show L < S^j^ > for the coupling 
constants J = 0.3, 0.4, 0.5, 0.75, 1, 1.5, 1.8, 2, 2.5, 3, 3.5, 4, 10 as a function of L/^k, beginning 
with L = 2. The data can be collapsed onto a single curve thus determining ^k{J) up to 
a multiplicative constant. If we fix one of the correlation lengths the rest of the correlation 
lengths are fixed by requiring that the scaling form be obeyed. An excellent data collapse 
is obtained. All the results in this figure are for the one impurity case in the ground-state 
subspace RPH = 1, = 1/2, with m = 128 states. In principle it is possible to obtain the 
Kondo length scale, C,k as a function of J from this scaling plot. However, since C,k varies 
quite rapidly with J one can essentially only obtain qualitative results at weak couplings 
where C,k is several thousand lattice spacings or at strong couplings where C,k is very small. 
However, for a fair range of intermediate couplings C,k can be extracted with a reasonable 
precision. Our results are summarized in Table |. We note that Fig. | clearly displays the 
complete crossover from weak to strong coupling. 

In the same manner we can look at L < Sj >, but instead of fixing j/L we can fix 
L/^k{J)- The scaling form should then be h{j /^xiJ)), where we have suppressed the 
dependence on L/^k{.J)- Using S,k{.J = 2.5) = 1.0, ^k^J = 3) = 0.6, such a data collapse is 
shown in Fig. |^ with L = 30, 20 for J = 2.5, 3, respectively. Thus, L/^k ~ 32 is kept fixed. 
Viewing ^k{J = 2.5) and ^k{J = 3) as fixed from the previous analysis this plot contains 
no free parameters. The collapse is excellent. The same rescaling can also be performed at 
other couplings. In Fig. |l^ we show results for J = 1.5 and J = 1.8 with C,k{J = 1-5) ~ 4.85 
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and ^k{J = 1-8) ~ 2.7, where we use L = 36, 20 for J = 1.5, 1.8 respectively. Again we find 
an excellent collapse of the data. The complete scaling shown in Figs. and ITO| is a highly 



non-trivial test of the scaling form Eq. (|8.1|) and the fact that the numerical results clearly 



follow the scaling form lends strong support to the existence of the Kondo length scale. 

B. Scaling of L < Sf^p > 

We now turn to a discussion of our results for the expectation value of the impurity spin. 



InFig.lyweshowL < 5f,^p > as a function of L for J = 0.5,1,1.5,1.8,1.9,2,2.5,3,3.5,4,10, 
beginning with L = 2. As is clearly evident L < S^^^^ > approaches a constant in the strong 
coupling limit. At weak coupling the behavior is consistent with L < S^^^ >~ L/2 for the 



values of L accessible. From the results in Fig. |TT] we can also check the result Eq. ( [7.141) . 



For J = 10.0 we find with the DMRG method that L < S^^^ >~ 0.0456 for large L. This 
can be compared to 0.0444 from Eq. ( |7.14| ), in very good agreement. 



C. 

In the previous section we obtained numerical results for C,k{J) by requiring our numerical 
results to scale. We can now try to fit these results to expressions for C,k obtained from 
renormalization group arguments which to first order gives Eq. ( p.lO|) . Is this possible ? 



The first point to realize is that we can only expect the weak coupling RG formula to work 
for a range of J such that 1 ^ ^ L. If C,k is too small then the coupling constant is too 
big so low order perturbation theory doesn't work. If C,k is too big then finite size effects 
will dominate. 

We proceed by obtaining a higher order expression for ^k{J) than the simple exponential 
relation Eq. ( |2.10| ). To go to one higher order, we first of all need to calculate the renor- 
malized coupling to 0[(Jp)^], obtained by reducing the effective bandwidth (in momentum 
space) from vr/a to 2/^0, where ^o, is some length scale much bigger than a, the lattice 
spacing. This involves the integral: 

29 



L ^ « ln|ta„(2&/a)l. (8.4) 

Thus: 

Ao = Jp + {Jpf ln[tan(2^o/a)] • (8.5) 
Here p is the density of states, 

p = l/2nt. (8.6) 

Next we need to integrate the /3-function, to third order. The /3-function isS: 

-dA/dlnA = - AV2. (8.7) 

We now integrate this equation, using the bare cut-off, Aq = l/C,o, the bare couphng Aq 
given above, the renormahzed couphng some number c of 0(1) and the new cut-off 1/C,k, 
the inverse correlation length. This gives 



l/Ao-l/c 



1 - (8.8) 



\ 1-2/c- 

This, together with the equation above for Aq determine vs. J to a better accuracy than 
the simple exponential form. We now have 2 free parameters; c which should be positive 
and 0(1) and which should be ^ 1 but <^ ^x- This form should be valid for the range 
of J where 1 <^ .^x <^ L. The corrections arising from the constant c is presumably only 
one of several equally important terms. We have included it here to improve the agreement 
with the numerical results. In the limit Aq Eq. ( p.8| ) reduces to: oc e^^^°/y/\o, in 



agreement with Ref. BO. 



A least square fit of the numerical results for C,k{J), in Table |, to the form Eq. ( pl8| 



is shown in Fig. |T^. The circles indicate the DMRG results while the solid line is Eq. ( ^.8|) 
with the fitted parameters = 1.76, c = 0.21. Taking into account the sizable uncertainty 
in C,K at weak coupling the fit works extraordinarily well. 
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IX. CONCLUSION 



We have shown that the local susceptibility can be brought into a unified scaling form 
linking a high temperature RRKY form and a low temperature local Fermi liquid form. 
In this picture the Knight shift is longer range at low temperatures where the screening 
cloud has formed, than at high temperatures where it has not. The experimentally observed 
behavior of the Knight shift0 can possibly be explained by a factorization of the scaling 
functions deep inside the screening cloud which is the only region the experiments probe. 

Non-trivial tests of the scaling form have been performed by numerically calculating the 
on-site magnetization at zero temperature for finite systems. In all cases scaling behavior 
consistent with the proposed form is observed. Both < S'£/2 > as a function of L, and 
< Sj > as a function of j for fixed C,k/ J clearly show scaling. Numerical estimates of Ck{J) 
has been extracted in good agreement with an estimate oiC,K{J) from renormalization group 
calculations. 
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TABLES 



TABLE I. The values for the Kondo length, ^^k-, used in Fig. |^ and shown in Fig. |T2| (circles). 
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W) 


10.0 


0.031 
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0.27 


3.0 


0.6 
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2.0 
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1.8 


2.7 


1.5 


4.85 


1.0 


23 


0.75 


100 
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1200 


0.4 


7000 


0.3 


100000 



35 



FIGURES 



CO 

V 



0.004 



0.002 - 



0.000 - 



-0.002 - 



-0.006 



1 1 1 1 1 1 1 T" 

- ® 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 


■ l\ !\ f ^ 




■ j 1/ y ^ 

1 ® 


y V V ••«>-^-^-^<»-^ «>- 


■0 





10 



15 
j 



20 



25 



30 



FIG. 1. The z-component of the electron spin as a function of site index j, for a chain of length 
L = 30. At the left end of the chain is an 5 = 1/2 impurity spin. The Kondo coupling is J = 0.05. 
The crosses indicate the first order perturbation result, Eq. (|6.5| ), the circles denote DMRG results 
for the state Sj< = 1/2, RPH = 1. m = 128 states were kept. The dotted line is the asymptotic 
result Eq. (U). 
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FIG. 2. The z-component of the electron spin as a function of site index j, for a chain of length 
L = 60. Impurity spins are present at both ends of the chain. The Kondo coupling is J = 0.05. 
The crosses indicate the first order perturbation result, Eq. ( |6.7| ), the circles denote DMRG results 
for the state = 1,RPH = 1, P = 1. 
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FIG. 3. The z-component of the impurity spin as a function of chain length, L. One impurity 
is present at the left end of the chain. The Kondo coupling is J = 0.05. The crosses indi- 
cate the first order perturbation result, Eq. ( |6.8| ), the circles denote DMRG results for the state 

= 1/2, RPH = 1. 
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FIG. 4. The correlation function < S^^^Sj > as a function of site index j (circles), for a chain 
of length L = 50. At both ends of the chain are S = 1/2 impurity spins. Only half the chain is 
shown. The Kondo coupling is J = 0.050, m = 150 states are kept with RPH = 1,P = 1, and 

= 1. The circles denote the DMRG results. The crosses indicate the perturbation result, 1/2 



times Eq. (6.5). 
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FIG. 5. The 8 states for the impurity site. 
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FIG. 6. The z-component of the electron spin as a function of site index j (circles), for a chain 
of length L = 24. At the left end of the chain is a 5" = 1/2 impurity spin. The Kondo coupling is 
J = 10, m = 128 states are kept with RPH = 1, = 1/2. The square denote the DMRG result 
for the impurity spin. The crosses indicate the perturbation result, Eqs. (^.6|), ( |7.11| ) and ( f7.16 ). 
The plus is the result for < 5imp > from Eq. ( 7.14| ). 



41 



A_ 

N 

CO 



V 



0.002 

0.001 

0.000 
-0.001 
-0.002 [ 
-0.003 - 
-0.004 \- 
-0.005 



I I I I 



I I I 



I I I I 



X 

O 



10 



15 



J 



FIG. 7. The correlation function < S^^^Sj > as a function of site index j (circles), for a chain 
of length L = 30. At both ends of the chain are S = 1/2 impurity spins. Only half the chain 
is shown. The Kondo coupling is J = 10, m = 150 states are kept with RPH = 1,P = — 1, 
and Sj- = 0. The circles denote the DMRG results. The crosses indicate the perturbation result, 
Eqs. (|7l9D and (|7^). 
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FIG. 8. Logaritmic plot of L < 5'£^2 > as a function of chain length L/(k for a range of 
different coupling constants. The initial point corresponds in all cases to L = 4. The solid lines 
are guides to the eye. The strong coupling limit corresponds to L < Sf^^^ >~ 1- 
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FIG. 9. L times the expectation value of the z-component of the electron spin, < Sj >, 
as a function of j/^K{J)- Two systems are shown: J = 2.5, = 1.0, L = 30 and 
J = 3.0, = 0.6, L = 18 Thus in both cases we have L/^^k 30. Clearly the data collapses onto 
a universal curve. 
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FIG. 10. L times the expectation value of the z-component of the electron spin, < Sj >, 
as a function of j/^K^J)- Two systems are shown: J = 1.8, = 2.7, L = 20 and 
J = 1.5, = 4.85, L = 36 Thus in both cases we have L/£^k ~ 7.4. Clearly the data col- 
lapses onto a universal curve. 
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FIG. 11. The scaled z-component of the impurity spin, L < S?^^ >, as a function of chain 
length, L, for a range of Kondo couplings J. In all cases does the first point correspond to L = 2. 
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FIG. 12. The Kondo length, as a function of Kondo couphng, J. The circles denote the 
numerical DMRG results. The solid line indicates a least square fit of the results shown to the 
form Eq. ( |8.8p . The fitted parameters are = 1-76, c = 0.21. 
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